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c basis  of  the  detailed  I 


problem,  and  the  second  method  approximates  the  photon  transport  through 


a specialized  Monte  Carlo  scheme.  The  results 


les;  additionally,  region  of  interest  images 


reduces  the  capability  for  feature  detection  and  analysis.  To  overcome 


le  extrapolation  algorithm  that  uaea  the  growth  or 
decay  oC  transform  coefficients  to  generate  sub-pixel  information  has 

extrapolating  under  sampled  LMBR  and  medical  x-ray  images. 


Lateral  Migration 


Radiography  IIMBR1 was  developed  at  the  Univeraity  of  Florida.  One 
auch  system  has  been  applied  to  detect  plastic  landatines  buried  in  soil. 
As  Che  name  suggests,  this  modality  uses  the  lateral  transport  of 


e strength  and  image  acquisition 


medium  as  well  as  the  surface  i 


collisions,  multiple-collision 
electron  density  variations  whi 

and  quantitatively  the  reasons  for  the  image  contras! 


density  and  identify  such 


to  lateral  photon  migration.  On  the  basis  of  the  photon  transport 


simulation  algorithm  (1 


from  low-resolution  images  accurately  is  difficult.  To  facilitate  the 
feature  extraction  from  such  low-resolution  Images  ovetcomplete  frame 

le  development  of  a new  algorithm  for  image  extrapolation  based  on 


AoW  algorithm  J 
modalities 


le  next  resolution  scale  without  blurring  the  edges  and 
! with  outstanding  contrast  details.  Additionally#  this 
Iso  helps  in  multi-resolution  feature  enhancement.  The 
applied  to  both  1MBR  images  and  images  from  other 

ure  enhancement  algorithms  have  a wide  area  of  application  ranging 
i remote  sensing  to  digital  video  image  enlargement. 

The  various  applications  of  CBI  along  with  their  features  are 
ussed  in  Chapter  2.  Chapter  3 presents  the  theory  behind  the  Monte 
o technique  for  studying  photon  transport  and  discusses  the  general 


simulation  algorithm 


OVERVIEW  OF  SCATTER 


production.  Additionally,  there  are  some  minor  effects  such  as  Raleigh 
nuclear  resonance  scattering  and  photodisintegration  of  nuclei.  Quantum 


sum  of  the  cross  sections  for  the  three  different  processes, 
transferred  to  the  bound  electron  » 

interaction  probability  increases  with  atomic  number  <E)  of 
material  and  varies  approximately  as  2s.  This  interaction  i: 

probability  at  energies  above  its  threshold  at  1.02  MeV  and 


relativistic  billiard  ba 


1 + att  - A) 


The  COKSCAH*  application  uses  Comptoi 

information  for  an  object  ia  generated.  To  facilitate  good  atatiatica  in 
signals,  both  aperture  and  detectora  are  kept  near  the  aurface.  Figure 

images.  This  essentially  limits  the  depth-of-view  in  the  material  to  a 
photon  mean  free  path.  This  modality  can  be  applied  only  to  image  smooth 
surfaces  because  multiple-collisions  from  the  surface  would  destroy  the 


Figure  2-2  Scatter  imaging  applications,  (a)COMSCAN  (b)X-ray  scatter 


material  along  chord  length  z.  Ao  la  the  optical  path  length  for  primary 
scattered  photons  and  A.  is  the  optical  path  length  for  the  multiple- 


each  detector  signal,  the  density  profile  along  c 


Raster  Gap 


ic  Individual  photon  history  is  simulated  a; 


incident  angles  ate  irrelevant  for  calculating 
excessive  computer  time.  Additionally,  mat 


interaction  probabilities 


Che  new  position  is  determined  through  the  scheme  as  described  by  Figure 
Figure  3-4  provides  a method  For  identifying  the  interacting 

many  commercially  avaiiable  computer  programs  such  as  MCNP15  (Monte 

flexibility  in  geometry  description,  source  description  and  material 
definition.  Because  of  this  flexibility,  the  MCNP  program  has  been 
chosen  to  simulate  the  particle  transport  in  the  U4BR  system. 

migration  backscatter  imaging  systems,  the  validity  of  MCNP  results  for 


Verification  of  the  Suitability  of  tf 


adiation  transport  code  uses  a sequence  of 
from  a known  random  number  sequence  and  a 


here  S indicates  the  stride,  R is  the  random  I 
radiation  transport  problems. 

generation  from  a tungsten  target  and  the  photon  transport  problem 
problems  are  quite  complicated  and  involve  simulation  of  detailed 


additionally. 


billion  his 


altar  this  period.  Thus,  MCNP  results 


u* 

“ T 

*-  252“ 

n MCNP  electron- 


The  results  in  the  Table  3-2  show  that  the  MCHP  predictions  have 


Figure  3-6  Experimental 


rectified 


H=  I 


thus  always  fixed  relative  to  each  other  instead  of  experiencing  a 


result  used  for  MCNP  verification  is  an  image  of  a 12"  diameter  plastic 


The  simulation  for  1"  per  foot  soil  slope 

for  the  MCNP  collimated  detect. 

for  a 1"  per  foot  slope  in  the  raster  direction.  The  HCNP  simulation  for 

pixels  and  the  simulation  was  repeated  for  the  next  16  pixels  to 

and  76  for  the  collimated  and  the  uncollimated  detectors,  respectively. 


:ks  plotted  by  the  original  SABRINA  represented 
tr  order  tracks).  The  original  program  was 
,s  modification  facilitated  the  study  of 


slightly  modified  to  al. 


jrce  positioned  o' 


MCNP  simulation  uses  200,000  source  particles  and  a 175  kVp  spectrum. 


Additionally,  most  of  the  scatter 


nd  1168  histories/6099  t 


first-collision  components.  This  imaging  modality  is  referred  to  as 
Lateral  Migration  Backacatter  Radiography  l LMBRI  because  the  image 
contrast  in  the  collimated  detectors  is  generated  by  a higher  lateral 


Verification  o: 


ABRINA  Results  with  Analytical  Calculations 


a(K-N)  cross-section  relations  are  applied, 
scattering  probabilities,  the  average  values 
angles  for  scattered  photons  are  calculated. 


urce  energy  that  yields  an  approximately  equi-distribution  o] 


photons . 


““  g.ss‘s  ;s'Sa;Ts,-Lr^sss™  — 


deposition  in  the  detectors,  the  desigi 


SOURCE  SPECTRA  SIMULATION 


5.1  E.N(E) -dE-Const.  Z.(T-E).dE 

N (El dE  is  the  number  of  photons  generated  for  a finite  energy  interval 


target  and  assumed  a perpendicular  incidence  of  the  electron 
beam  on  the  target.  A detailed  physics  based  model  was  recently 
published  by  Tan  and  Heaton  ” and  their  method  models  both  an  angular 


Heaton,  respectively.  The  characteristic  x-rays  Chat  are  expected  to 


MM 


5.5  p.  CM  = H’  - T‘ 

5 6 rain  ■ .xp^^SsLzJ^l Co3t>i j 


-ng  energies  and  fractional  yields  considered 


t,  -l)"/(Ojp(|).-"‘ 


purity  germanium 


fl-Xb/Xcl, 


175kVp  and  200  kVp  spectra,  respectively.  1 
with  measured  results  is  better  than  for  the  XSPEC  results  because  MCNP 


i MIGRATION  BACKSCATTER  {I 


The  MCNP  Monte  Carlo  analysis  in  Chapter  4 has  shown  that  the 
images  formed  in  the  collimated  detectors  are  predominantly  due  to 


A two-collision  approximation  was  developed  by  Hatanabe  et 


fj(*y!*o.y)  - c|i-»«»i(-<^(ai  + " + D){«pC3(UnyiMi(],y0))-i)j+n 


x-cay  source  projection  to  the  pixel  projection 


F,(x,  y,  x„,  y«)  = c,(l  + exp[-c,(h  + hi  + s,)  ] lexp|C,l,)  - 11)  + 17 


=(h+H/2)/cos (9j) 
-1/cos  (6i)  and 
■s/co5(8i) . 


6.6  s,(x,y)»s,(x,y-2yo)  • 


8 (D,  - D,)‘  - [s’  + 2as,l*,  y)  - S,(*  yl )]  - 45,  (x,  y|s,(x,  y)  + [s,(x,  y)  + s,lx,  y)]' 
ere  £ = 0,1*,  y)  - fl,lx,y). 


9 (D,  -D, 


■ [=,(*.  y)  + Sjlx,  yl]’  - 4 5,0,  yBjIn,  yl 


Recasting  tl 

6.10  [3,1*,  yl  + S,lx,  y)]’  = [d,(x,  yl  - D,(x.  y)]’  + [4S,(x,  yls,!*,  y)]1" 

where  Che  '5'  on  the  FHS  is  approximated  by  the  two-collision  model  ai 


dified  two-collision  model  ar< 


Specialized  Monte  Carlo  Hodel,  Photon  Lateral  Migration  Simulator ( PLUS) 


MCHP  terminates  a particle  history  according  to  particle  weight  and 


collision,  direction  of  travel  or  migration  length.  The  MCHP  allows 


relevant  to  the  problem.  This  automatically  forces  the  use  of 


4.  Additionally,  the  detailed  physics  modeling  ar 


whether  the  particle  is  in  mine  or  soil,  and  supplies  a 
flag  for  variable  'mat'  indicating  the  type  of  material 


- Implements  implicit  capture  and  applies  the  weight 
biasing  factor  'f'  and  calculates  the  weight  *wei'  i 


it  scattering  c 


le  mean  and  standard  deviation  fi 


e the  surface  is  tailied 


simulation  used  90,000  particles  per  pixel  to  generate  a 33  X 33  image. 


signals  are  predicted  with  an  uncertainty  of  10%.  The  trend  lines  fi 
Figure  6-6  and  the  trends  of  measured  and  MCNP  results  presented  in 
Figure  3-14  have  similar  behavior,  even  though  the  uncertainties  f< 
the  collimated  and  uncollimated  detector  signals  in  the  latter  case 

generated  by  the  PIUS  model  for  LMB  photon  transport  are  a good 

The  PLMS  collimated  detector  simulations  took  65  minutes  on  a 
Pentium  90  MHz  personal  computer  for  generating  each  33X33  image  Mi' 


o-parallel  option  the  simulation  t: 


<0  photons  per  pixel  at  relative 

computer  with  very  low  mine  image  contrast  and  requires  at  least  200 
particles  to  generate  the  same  image  contrast  as  PLMS.  PLMS  due  to  it 

second-collision  photons  and  hence  produces  a high  contrast  LMBR  min 
inferred  that  PLMS  also  aids  in  speedy  simulation  of  LMB  systems. 


Figure  6-5  PLHS  generated  images  for  the  front  collimated  detector,  (a) 
the  12'  plastic  landmine  at  (23,23),  (b)  12"  metal  »ine  dob'^LHs'  *"d 
90,000  particles/pixel  and  the  pixel  intensities  haue  uncertainty  of 


MU 


■ (inches) 


artifacts.  Low-resolution  images  can  broadly  be  classified  as  under 
sampled  images.  Wavelet  based  multiresolution  analyses  are  less 

reduces  the  capability  for  feature  detection  and  analysis."  To 
frame  representation.  An  algorithm  tt 


implemented.  An  imag 


1 of  wavelets (AoW)  was 


Imaging  modalities  such  aa  lateral  migration  b. 
radiography  (LMBR)  produce  1< 
under  sampling  in  such  imaging  m 

requirements  o£  the  imaging  modality,  and/or  poor  lightii 

interest  selected  from  a high  resolution  image  o: 
sampled  image  when  magnified  via  traditional  digital 

resolved  hence  making  it  difficult  to  establish  a correspondence 


te  lies  in  the  local  variation  of 
mage  the  discontinuity  in  pixel 

epresentations  of  image  information  under  various 
:h  orthonormal  wavelets  (that  are  periodic 

wavelet  expansion  of  any  vector  without  sub-sampling  would  also  contain 
circularly  shifted  wavelet  coefficients  of  the  original  signal.  This 
periodic  behavior  of  the  wavelet  coefficients  is  due  to  the  asymmetric 
nature  of  standard  wavelets  used  in  image  c< 


some  loss  of  information.  The  projection  of  a function  using  any  ba 

signal  is  shift  invariant,  it  facilitates  the  study  of  growth  or  d 
of  the  image  projections  across  resolution  scales.  By  evaluating  th 


polation,  such  as 
interpolation,1'  pixel 


7-j  js - £.v* m i°> 


VJs  are 


The  two-dimensional  filters  o; 


■ ■EE 


Figure  7-S^The  log-log  plot  o£  filter  coefficients  for  32  pixels  verses 


A - ESTIMATED  L7  COEFFICIENT  B-  ACTUAL  L7  COEFFICIENT 


donates  the  scale  size  of  the  detail  function. 


detail  coefficients . 


Obtain  the  image  projections  Ln  in  the  Fourier  domain  wit, 

AoH  shell  coefficients  in  the  spatial  domain. 

For  every  pixel  position , use  the  pixel  intensities  for  a. 


expanded  image 


original  image. 


implemented  in  Matlab.  30  To  understand  the  effect  of  noise  on  the 
multiresolution  analysis  and  the  image  sub-pixel  information  generation, 


Comparison  of  AoW  approximation  with  other  interpolation  methods 


AoW  projection  algorithm  is  compared 
Figure  7-9,  and  the  mean  squared  erro 


extrapolation  approximates  an  image 


the  original  henna  image  in 


transform  based  image 


techniques  such  as  pixel  replication,  line 
interpolation,  pixel  regularization33  and  I 


Ith  other  established 


led  portion  <256  by  256)  of  the  henna  image.  The  interpolated 


original  image 


extrapolated  image 


Figure  7-10  Comparison  of  image  interpolation  algorithms (a)  the  original 
(b)  pixel  replication  (c>  bi-linear  interpolation  <d)  bi-cubic 
interpolation  (e)  pixel  regularization  (f>  b-spline  interpolation  (g) 


the  cubic  - interpolation  gives  smooth  results,  whereas  B-spline  and 


interpolated  image  with  respect  to  sharpness. 


are  many  modalities  th 
:h  modality,  LMER,  is  VI 


y promising.  However,  the 
under  sampled  and  detecting 


The  zero-crossing  representations  of 
landmine  field  image  is  shown  in  Figure  7-3 

show  the  details  of  a circular  step  patten 

the  rectangular  object  is  indeed  a mine.  Additionally  the  ring  patterns 
detection  ability  greatly  assists  in  the  identification  of  these  LMBR 


y difficult  tc 


r coefficients  o: 
er  applying  the  i 


5 

20 

30 


*9 


typical  cylindrical 


of  an  image  computed  from  autocorrelation 


wavelet  transform  coefficients 
shell  of  an  wavelet  can 


the  statistical  basis  of  the  trend  of  AoW  coefficients  at  each  pixel, 
sub-pixel  information  is  approximated  at  a higher  resolution.  The 
quality  of  an  image  expanded  through  the  transform-based  AoW  shell 
algorithm  is  comparable  to  or  exceeds  the  quality  of  an  image  expanded 

mathematical  framework  for  carrying  out  multiresolution  analysis  and 


Additionally,  l 


capability  of  the 
nhancing  LMBR  images.  Though  the  LMBR  system  studied 


. The  Monte  Carlo  method  of  simulating  photon  transport,  specifically 
photon  transport.  The  random  number  generator  used  in  MCNP  is  found 

photon  transport  problems  is  as  great  as  50%.  This  difference  is  due 
to  the  idealization  of  a complex  electron  source  in  the  MCNP 

ns.  Additionally,  MCNP  does  not  consider  long  range  e-e 

IR  based  landmine  imaging  system  has  been  analyzed  using 

electron  density  (low  z and  low  density)  allow  higher  lateral 


of  generating  two  images  one  containing  predominately  surface 

density  materials.  This  improves  the  ability  to  distinguish  plastic 
explosive  from  soil,  wood,  and  vegetation.  The  photon  lateral 
migration  length  is  different  for  each  material.  The  difference  in 
lateral  migration  length  leads  to  unique  signatures  for  different 

signatures  can  be  used  to  approximately  determine  the  material  of  the 
imaged  object.  This  feature  discrimination  capability  of  LMBR  is 


differentiating  plastic 


The  one-dimensional  x-ray  spectrum  prediction 
developed  and  its  results  compared 


simplified  collision- 


establishing  the  theory  dictating  a new  imaging  modality,  LMBR  has  been 
developed  along  with  a new  transform  based  approach  to  the  two- 


dimensional  interpolation  problem, 
developed  here,  many  applications 
extrapolation  can  be  developed.  ’ 


i.  The  two-collision  model  developed  here  is  very  specific  to 

cylindricel  plastic  objects;  its  applicability  to  other  scattering 
materials  has  to  be  verified.  Additionally,  a two-collision  kernel 
along  with  image  subtraction  can  be  used  for  correcting  the  multiple- 
scatter  photon  effects  in  other  scatter  radiography  techniques. 

>.  PI-MS  unlike  MCHP  includes  the  collision  termination  feature.  This 
particle  termination  efficiently  eliminates  the  influence  of  f: 


An  artificial  intelligence (All  algorithm  to  discriminate  different 
objects  in  LMBR  images  can  be  developed.  The  AI  algorithm  can  use 


objects  and  can  be  trained  through  a set  of  MCNP  c 
chapter  7 is  in  dyadic  scales  and  con  be  extended  to  integer  sc 


TRANSPORT  EQUATION 


APPENDIX  A 

THE  MONTE  CARLO  SOLUTION  TO 

can  consider  either  the  integral  or  integro-differential  £• 
or  neither.  The  selection  o£  the  form  of  the  BTE  is  entirely  d« 
on  the  problem.  The  integral  form  of  the  BTE  will  be  considered  here 
because  the  MC  solution  to  the  BTE  in  this  form  is  straight  forward 
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Continuing  the  same  process  for  the  phase  volume  dR.  about  R.  and 


Equation  A. 7 has  a drawback  in  that  it  breaks 
acatterinq  cross-section  for  a qiven  material 


assumed  at  that  point  and  the  photon  weight  W is  suitably  adjusted 


e integral 


CALCULATION  OF  AVERAGE  ENERGY  AND  ANGLES  FOR  COMPTON  SCATTER 
The  fallowing  equations  are  used  to  calculate  the  average  scattering 


( 2.  *.  sin  $.  f.  dtf  . 

forward  scattering  cross  section  is  obtained  through 
! Z a.  sin  0.  t.  dO  . 

backscattering  crosa  section  is  obtained  from 


= J Zc  sin  6.  cos  0.  f.  d 0 ar 


B6  cos  0_  = — j 2x  sin  ft  cos  ft  t.  dO  . 
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